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Using the analogy with inelastic granular gasses we introduce a model for wealth exchange in 
society. The dynamics is governed by a kinetic equation, which allows for self-similar solutions. The 
scaling function has a power-law tail, the exponent being given by a transcendental equation. In 
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I. INTRODUCTION 

The distribution of wealth among individuals within a 
society was one of the first "natural laws" of economics 
PJ. Indeed, its study was motivated by the desire to 
bring the accuracy attributed to natural sciences, namely 
physics, to economic sciences. The celebrated Pareto law 
states that the higher end of the wealth distribution fol- 
lows a power-law P(W) ~ W~ 1 ~ a with exponent a ro- 
bust in time. 

The validity of the Pareto law was questioned and re- 
examined many times but the core message, stating that 
the tail of the distribution is a power law remains in 
force. There are recent investigations, e. g. Q,H,B,|^, 
giving reasonable empirical evidence for it. In fact, it 
is not so much the functional form itself but its spatial 
and temporal stability that is intriguing. Indeed, while 
the value of the exponent a may slightly vary from one 
society to another, the very fact of the power-law O.tail 
in the distribution is valid almost everywhere. Recent 
investigations suggest that the range of validity of the 
Pareto law may extend as far in the past as to the ancient 
Egypt of the Pharaohs . 

The universality of the power-law tail is surely a phe- 
nomenon asking for explanation. Recently, there was a 
lot of effort establishing finally the multiplicative random 
processes repelled from zero as a mathematical source of 
the power-law distributions 0, IE 0, E3, EH E3, E3, IS EE 
EE IlE EE EE E3|. Alternatively, the killed multiplica- 
tive processes as sources of power-laws were studied in 
U However, there are plenty of possible ways how the 
multiplicative random processes of this type come onto 
scene. One of the most studied implementations were 
the generalized Lotka-Volterra equations 0, [ill E3, 
and the analogy with directed polymers in random me- 
dia 0, I22I [23 ■ Both of these schemes are formalized 
by a kinetic equation describing the exchange of wealth 
between agents and global redistribution of wealth which 
plays the role of repelling from zero. Related approaches 
were subsequen tly pur sued by a number of studies and 

simulations 0, M, M ^m^MM,M,M,M,M, 
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More recently, empirical studies of the lower end of 
the wealth axis showed that the distribution of wealth is 
rather exponential than power-law, while the high- wealth 
tail still remains power-law 0, EI D3 • This finding was 
interpreted as a result of a conservation law for total 
wealth, leading to the robust Boltzmann-like exponential 
distribution, whatever the random wealth exchange be, 
in full analogy with the energy distribution in a gas of 
elastically scattering molecules. 

This, together with older studies within the same spirit 
[iflj . lead to the view of economic activity as a scatter- 
ing process of agents, analo gou s to inelastically scattering 
particles p^ . I30L l3ll EI lIolTiri lirj . Indeed, the inelastic- 
ity is indispensable to explain the power-law tail and it is 
also reasonable to suppose that the total wealth increases 
on average. 

The numerical simulations performed to date confirm 
the emergence of power-law tail in agent-scattering pro- 
cesses with great reliability. However, analytic insight is 
lacking in most of the studies available today. The main 
concern of our work is to fill this gap, providing ana- 
lytical results at least for a simplified model of wealth 
exchange. To comply with the task we will be guided by 
existing analytical approaches for models of inelastically 
scattering particles. 

Inelastic scattering of particles was studied thoroughly 
in the context of granular materials ^E|- The simplest 
one of the models used is the Maxwell model , w hose in- 
elastic variant was investigated in detail [43, [5(J, OHl 152 . 

HHmmi^ilililiil. More realistic mod " 

els of granular gasses were also introduced [iE E3| but 
their full account goes beyond the topic of this work. 
The most important conclusion of these studies is that a 
self-similar solution of the kinetic equations exist, which 
is not stationary in time, but assumes time-independent 
form after proper rescaling of the energy. The tail of the 
scaling function becomes power-law under certain condi- 
tion. 

The formalism developed for granular gases can be 
readily adapted for binary wealth exchange of agents. In- 
deed, within the mean-field version of the Maxwell model 
the particles scatter randomly one with another irrespec- 
tively of their positions. This corresponds to randomly 
picking pairs of agents for interaction, with no care of the 
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(possibly complex) structure of their relationships. In re- 
ality the economic activity goes along links in a complex 
social network Indeed, recently there were inves- 

tigations of the role of network topology in wealth distri- 
bution [34l |(5(| . We may consider the present model as 
an approximation of that network by a complete graph. 

The main difference from the mean-field Maxwell 
model is that the energy of the granular gas decreases by 
dissipation, while the average total wealth of the agents 
increases due to the economic activity. The sign of the 
non-conservation is therefore opposite in the two cases. 
While the form of the equations may remain the same, 
the solution cannot be directly continued from one do- 
main to another. Therefore, while the case of dissipation 
is relatively well understood, new approaches are needed 
in the case of production. That is the aim of the present 
work. 



II. INTERACTING AGENTS AS SCATTERING 
PARTICLES 

A. Description of the process 

Imagine a society of N agents, each of which possess 
certain wealth m, z = 1,2, N. Time-to time the agents 
interact in essentially instantaneous "collision" events, 
when certain fraction of the wealth can be exchanged. 
Moreover, we suppose the system is open and the inter- 
action can catalyze an increase of the total wealth of the 
two interacting agents. Indeed, the source of the human 
wealth lies beyond our society and the ultimate cause is 
the energy poured to the Earth from the Sun. Nonethe- 
less, the external energy is utilized only through a human 
activity and we simplify the problem by assuming that 
the net increase of wealth happens at the very moments 
of agents' interaction. 

We also assume that only pairwise interaction occurs. 
This may be a very crude assumption, as corporate de- 
cisions affect many agents simultaneously. However, we 
expect the presence of multilateral interactions does not 
affect the essential mechanisms in work here. 

The dynamics of our model is described as follows. In 
each time step t a pair of agents is chosen randomly. 
They interact and exchange wealth according to the sym- 
metric rule 
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All other agents leave their wealth unchanged, Vk(t + 
1) = Vk{t) for all k different from both i and j. The 
parameter f3 £ (0, 1) quantifies the wealth exchanged, 
while e > measures the flow of wealth from the outside. 
The process is sketched schematically in Fig. [TJ 

This rule is similar to those studied in ^3|> [53, Ei| and 
simulated numerically in [2^ . I3TL lil . |4^| but we consider 
it slightly more realistic as it treats the agents in a priori 
symmetric manner. It also embraces various sources of 




FIG. 1: Schematic picture of the scattering process, where 
the wealth is exchanged and produced. 



wealth non-conservation within a single effective parame- 
ter e. In fact, also the formulation based on the similarity 
with the problem of directed polymers 0, |23| can be re- 
duced to a rule of the form similar to Q). Therefore, 
we are studying a representative of a whole class of re- 
lated models and we expect the analytical results we will 
present have rather broad relevance. 



B. Kinetic equation 

The equation describes a matrix multiplicative 
stochastic process of vector variable v(t) in discrete time 
t. Processes of this type are thoroughly studied e. g. in 
the context of granular gasses. Indeed, if the variables u, 
are interpreted as energies corresponding to i-th granular 
particle, we can map the process to the mean-field limit 
of the Maxwell model of inelastic particles. However, the 
energy dissipation conventionally quantified by the resti- 
tution coefficient implies now the negative value e < 0, 
contrary to our assumption e > 0. We will see later that 
this apparently small variation makes big difference in 
the analytical treatment of the process. 

The full information about the process in time t is 
contained in the iV-particle joint probability distribu- 
tion Pjv(t; V\, i>2, Vn)- However, we can write a kinetic 
equation involving only one- and two-particle distribu- 
tion functions 



P 1 (t+l;v)-P 1 (t;v) = — 



-Pi (*;«)- 



+ / P 2 {t; Vi,Vj) S((l -{3 + e)Vi + pvj - v) dvid Vj 



(2) 



which may be continued to give eventually an infinite 
hierarchy of equations of BBGKY type. As a standard 
approximation we use the factorization 



P2{t;v u v j ) = Pi(i;u i )P 1 (t;v i ) 



(3) 



which breaks the hierarchy on the lowest level, neglect- 
ing the correlations between the wealth of the agents, 
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induced by the scattering. In fact, this approximation 
becomes exact for N — ► oo. Therefore, in thermody- 
namic limit the one-particle distribution function bears 
all information. 

Rescaling the time as t = 2t/N in the thermodynamic 
limit TV — * oo, we obtain for the one-particle distribu- 
tion function P(t;v) = Pi(t,v) a Boltzmann-like kinetic 
equation 



dP(v) 

dr 



P(v) 



P(vi)P(vj) <5((1 -0 + e)vi + /3v 3 - v) dv z dv 3 



(4) 



which describes exactly the process Q in the limit N — > 
oo. This equation has the same form as the mean-field 
version for the well-studied Maxwell model of inelasti- 
cally scattering particles [5-iL l5rl l57j . The main differ- 
ence consists in the fact that here the wealth increases, 
while in inelastic gas the energy decreases. This seem- 
ingly little difference has, however, deep consequences for 
the solution of Eq. (@J. 

Note also that within the framework of Maxwell model 
the distributions are expressed in terms of velocities, 
while our dynamical variables correspond rather to en- 
ergies of the particles. 



III. SOLUTION OF THE KINETIC EQUATION 

A. Self-similar solutions 

Note first that the average wealth v = J v P(v) dv in 
the process described by the kinetic equation Q grows 
exponentially 



v(t) = v(0)e e 



(5) 



and therefore Eq. has no stationary solution. How- 
ever, we may look for a quasi-stationary self-similar so- 
lution in the form usei mil 



P(t;v) 



v(t) v{t) 



w e 



(6) 



•dw 



Using the Laplace transform = J °° $( 

we can write a non-local differential equation for the scal- 
ing function in the form 

ex&(x) + $(x) = $((1 -0 + e)x) 3>(0x) (7) 

A hint about possible solutions can be obtained from 
a special exactly solvable case e = — 2^/0 +2/3. It 
can be easily verified [Hi! that the function $1(2;) = 
(1 + \/2x)c~ v ^ is a solution of J7J). Inverting the Laplace 
transform we obtain the corresponding wealth distribu- 
tion <&\{w) = —L= w~ 5 / 2 exp(— -M which has similar 

form as obtained in previous studies UJlml^a. How- 
ever, in this case the value of e is negative, which contra- 
dicts our assumption of wealth increase, while for e > 



the above idea leading to the function $i(x) does not 
work. Therefore, we must look for alternative ways. The 
leading idea of our approach is that equation J7J is nearly 
local for small values of e and 0. Therefore, we will ex- 
pand the factors on the RHS of Eq. (0 in Taylor series 
in e and and perform the limit e,0 — > 0. As the pa- 
rameters e and quantify the amount of wealth increase 
and exchange in single trade event, we interpret the lat- 
ter limit as the limit of continuous trading. In fact, such 
limit should involve also a rescaling of time t, but because 
we are interested only in stationary regime, the explicit 
time dependence does not enter our considerations. 

It should be also stressed that an important feature 
can be inferred from the observation that the system be- 
haves differently for positive and negative e. Indeed, it 
suggests a singularity at the point of precise conservation 
of wealth, e = 0. 
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FIG. 2: Solution of the equation f (j3, e, a) = (1 + e - j3) a + 
(3 a - 1 - ea = for e = 0.1 and = 0.0025 (full line) and 
= 0.004 (dashed line). 



B. Power- law tails 

The main concern in empirical studies of wealth distri- 
bution is about the shape of tails, which assumes power- 
law form. The behavior of the distribution (fr(w) for 
w — ► 00 can be deduced from the singularity of the 
Laplace transform <t(x) at x — > 0. Therefore, we assume 
the following behavior 0, 



$(z) = 1 -x + A 



.. for x — ► 



(8) 



where a € (1, 2). This type of singularity results in the 
power-law tail as $>(w) ~ for w — > 00. Insertion 

of (JSJ) into J2J leads to a transcendental equation for the 
exponent a 



{l + e-0) a +0 a -l-ea = Q 



(9) 
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the solution of which is illustrated in Fig. |2J Obviously, 
there is always a trivial solution a = 1. The power-law 
tail is due to another, non-trivial solution, which falls 
into the desired interval (1, 2) only for certain values of 
the parameters (3 and e. We can see the allowed region 
in Fig. solution in the range a <E (1,2) exists within 
the shaded region. We can also see that fixed value of 
a defines a line in the (3-e plane. We can approach the 
limit e — > 0, (3 — > while keeping a constant. This is 
to be interpreted as continuous trading, as the amount 
of wealth exchange and increase in a single trading step 
is infinitesimally small. Making this, the non-local terms 
in Eq. J7|) become local and we can expect to obtain an 
ordinary differential equation, soluble by standard meth- 
ods. 
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FIG. 3: Solution in the range ot G (1,2) exists within the 
shaded region, dashed line corresponds to a = 2, dash-dotted 
line corresponds to a = 1 and full line to the solution a, = §. 



C. Continuous trading limit 

Indeed, expanding @ we obtain the following formula 
relating (3 and e for fixed a in the limit of continuous 
trading (3 — > 0, e — > 0: 



= 



a — 1 



+ 0{e A ) + O{e zo ) 



(10) 



The leading correction term to depends on the value 
of a; for 1 < a < 3/2 it is of order 0(e 2a ), for 3/2 < 
a < 2 it is of order 0(e 3 ), while in the special point 
q = 3/2 we should include both correction terms, as 
they are of the same order 0(e 3 ). Systematic expansion 
in e is developed in Appendix 1X1 

Taking the same limit with fixed a in Eq. we ob- 
tain, using 11U|) . the following equation 



1 



x$"(x) 



1 



. 



(11) 



Of the two independent solutions of (|11|) only one has 
correct asymptotics $(x) — > for x — > +oo. It can be 
expressed using modified Besscl function 



$(x) = C x a/2 K a (2 V^TVx) 



(12) 



where the constant C is fixed by the normalization 
$(0) = 1. Inverting the Laplace transform we finally 
obtain the wealth distribution 



$(w) = CuT 



a - 1 
exp( ) 



(13) 



with C= (a - l) a /T(a). 

We can see that the distribution obtained exhibits the 
desired power-law behavior for large wealth. Moreover, 
it has a maximum at a finite value of w = w max = 
(a — l)/(a+ 1) and depression for low wealth values. 
The size of the depletion is determined by the exponen- 
tial term in (|13|l . i. e. by the same value of a which de- 
termines the power in the power-law. This corresponds 
to the idea presented e. g. in Rcf. ^l| stating that it 
is the value of the lower bound for the allowed wealth 
which determines the value of the exponent. Here, how- 
ever, this result comes purely formally as a result of the 
analytic computation. In our approach it is the inter- 
play between wealth increase (parameter e) and wealth 
exchange (parameter (3) that dictates the value of the 
exponent a. 



D. Corrections for finite trading in one step 

Expanding the equation Q in powers of e and (3 it 
is possible to include systematic corrections to equation 
(|ll|l and therefore corrections to wealth distribution l|13|l 
for finite amount of wealth increase and exchange in sin- 
gle trading step. Details of the calculations are given in 
Appendix^ here we only summarize the results. 

The expansion (|l(Jf) of the parameter (3 in powers of e 
can be continued as 



= 



a-1 



1 fa - 1 
a 



e 2a - 



(q-l)(2a-l) c 3 (14) 



+ 0(e 4 ) + 0(e 4Q - 2 ) 



Correspondingly, the wealth distribution, expanded in 
powers of e is 



$(w) 



(a - 1) Q 



-l-a 



exp 



l-a 



T(a) \ w 

a — 1 (la a — 1 

1 - 



2 

a 



3 

a — 1 



^io e 



lnui 



1 

w 



.2(a-l) 



+ 0(e 4 ) + 0(e 4a ~ 2 ) 



(15) 
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where the constants isqi and viq are given in Appendix^ 
We show in Fig. 0|the wealth distribution according to 
(115(1 for a = 1.7 and several positive values of e, namely 
for e = 0.03, 0.1, and 0.3. We can see that the distribu- 
tion is affected mainly at small values of wealth, shifting 
the maximum toward smaller w when e increases. On the 
contrary, the tail of the distribution is nearly unaffected, 
showing universal and robust power-law behavior. 
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FIG. 4: Wealth distribution according to Eq. I|15|l for e — > 
(full line), e — 0.03 (dashed line), e = 0.1 (dash-dotted line), 
and e = 0.3 (dotted line). 

Let us stress again that the solution known for e < 
cannot be properly continued to the region of e > 0, due 
to the presence of singularity at e = 0. The singularity 
can be seeen e. g. in the behavior of the solution of 
Eq. @, as shown in Fig. [SJ However, for a = 3/2 the 
formula (|13fl describes the solution of J7J) on both limits 
e — > + and e — > 0~. This implies that the singularity 
is rather weak, because the solution of Eq. Q) is contin- 
uous in e, and only the derivative with respect of e has 
a jump at e = 0. One may speculate about the fate of 
the singularity if we allowed e and (3 not fixed param- 
eters but random processes themselves. Most probably 
the singularity would vanish but final answer is left for 
future work. 



IV. CONCLUSIONS 

We formulated a model of wealth production and ex- 
change, where agents randomly interact pairwise. Using 
the analogy with the mean-field version of the Maxwell 
model for inelastic scattering of granular particles we ob- 
tain analytical results for the wealth distribution. 

The dynamics of the model is governed by a kinetic 
equation for one-particle distribution function. We look 
for self-similar scaling solutions, corresponding to redefin- 
ing the unit of wealth after each wealth increase. The 
form of these solutions is given by a non-local differential 
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FIG. 5: Solution of equation © for a = 3/2 in the ranges 
e > (full line) and e < (dashed line). Note the singularity 
at e — which means that we must skip from one of the three 
solutions of @ to another one. 

equation, exactly soluble only in the practically irrele- 
vant case of net wealth decrease. Therefore we turned to 
approximation schemes. 

First, we looked at the behavior for large wealth. The 
tail of the wealth distribution has a power-law form, and 
its exponent a is determined by the interplay between 
the intensity of the wealth exchange and the amount of 
wealth produced. The form line in the /3-e plane with 
fixed a is found, depending quadratically on e for e — > 
0. The physically allowed values a G (1,2) determine a 
horn-shaped region in the (3-e plane. 

The second approximation consisted in taking the limit 
of continuous trading, meaning small wealth production 
and small exchange within a single trading operation, 
while keeping the exponent a constant. Here we ob- 
tained closed formula for the entire wealth distribution, 
which has power-law tail as expected and a maximum 
at certain (low) wealth value. The form of the wealth 
distribution corresponds to those found in previous stud- 
ies ^| 0, ■ It is interesting to note that this general 
form has one-to-one correspondence between the position 
u; max of the maximum of the distribution and the value of 
the exponent. There are few agents having wealth below 
■Wmax- This suggests that the intuition formalized c. g. in 
[Il|, [lJl, that the exponent is "tuned" by the low- wealth 
behavior of the distribution, may be in work quite gener- 
ally. Here, the free parameters are apparently the wealth 
production and exchange, but in reality these parameters 
may be themselves tuned by a mechanism which fixes the 
position of the maximum of the wealth distribution, i. c. 
the lowest wealth compatible with survival. 

However, there is still open question of the specific val- 
ues of the exponent, which arc quite robust in different 
societies. It seems, also on the basis of our results, that 
it cannot be explained by the bare mechanism of eco- 
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nomic exchange and some other ingredient, possibly of Agency of the Czech Republic, 
sociological origin, is required. 
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APPENDIX A: SYSTEMATIC EXPANSION FOR 
SMALL e AND j3. 

Let us start with the special value a = 3/2. Here, the 
equation Q has an explicit solution in the form 



J 



-3 17/3- 29/3 3 / 2 + 15/? 2 + 4/3 5/2 - 4/3 3 + ^^(3 - 2 y/ff) (3 (2 yffl + l) 3 - l) e 

v7J-3/? + 3/? 3 / 2 -/? 2 

I 



(Al) 



However, the non-local differential equation 0) still does 
not yield explicit solution. Inverting the expression (|A1|) 
we get the following series expansion 



a 1 2 1 3 1 4 7 5 113 

3 = -e 2 e 3 H C e 5 H < 

' 4 12 16 144 2592 



O 



(A2) 



For general value of a the variable f3 is expressed as a 
series in two small parameters e and r] = e 2 ( a_1 ), which 
coincide only if a = 3/2. Therefore, we can write 



E A- 



m+2(a— l)n 



(A3) 



and the various terms take variable precedence in the 
order of smallncss when e — > 0, depending on the value 
of a. For the first several coefficients we have 

Poo = ^ (A4) 



010 
A)l 



(a - l)(2a- 1) 
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(A5) 
(A6) 



Starting from the expansion l|A3|l we can convert the first 
order non-local differential equation J7J) for $(x) into 
infinite-order local differential equation for The 
price to pay for it is that the coefficients in the latter 
equation contain the moments fik = J $>(w)w k Aw of the 
solution itself. Indeed, we can write 



$((1 -P + e)x)= lim exp (e - (3) x— ) $(y) (A7) 
y ^x \ Ay; 



<&(f3x) = lim exp (3x— $(y) 
y^o V Ay; 



(A8) 



Therefore, we obtain a linear combination of terms of the 
following form 



l+n d m $(x) d»<i>(o) 



Ax r 



Ax r 



(A9) 



which, after inverse Laplace transform, give rise to terms 

Am+n 



(-1) 



Aw m+n 



[w m *(«;)] 



(A10) 



However, the first two moments are fixed by definition. 
Indeed, the normalization of the probability distribution 
fixes the zeroth moment and the fixed average wealth, 
imposed by the scaling condition JHJ fixes the first mo- 
ment, so that /l«o = Hi = 1. This consideration leads to 
the equations for lowest correction to the solution H13fl . 
which are free of unknown higher moments. 

Generally, the solution can be then expressed in the 
form of the series in powers of e and e 2 ( Q_1 ) 

oo 

$H=4>oM E <t>mnHe m+2 ( a - 1)n . (All) 

m,n— 

We assume (j)oo(w) = 1. The normalization must be in- 
dependent of e, which can be written as 



$o(w)4> mn (w)dw = S m0 S n0 . 
Therefore, the lowest term obeys the equation 



(A12) 



W XL) f \ a 



-w 



$o(w) = (A13) 



which has the following solution satisfying the normal- 
ization (| A 12(1 



(a 



T(a) 



1) Q -l-a I 1 
w exp 



(A14) 



Indeed, it coincides with the result of 1131) . 

The next two terms satisfy the following equations 



<t>'io( w ) 



a — 1 



1 



w 2 1 / a - 1 



(A15) 



(io-l) (A16) 
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which can be easily solved. We obtain where fy(x) = r'(x)/r(x) is the logarithmic derivative of 

the gamma function. 

, , -s a - 1 /2a a - 1 \ . . „ _ N 
<t>w(w) = — =- - i/io (A17) 



3 \ w w 2 
0oiH - ~{^ L ) (ln^ + ^-^i)(A18) 

and the constants i/qi, fio are fixed by the normalization 
condition (|A12(I . We find explicitly 

vio = a (A19) 

i/ i = In - * (a) + — — (A20) 

a — 1 
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